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A METHOD DISTINGUISHING LINE OF SIGHT (LOS) 
FROM NON-LINE OF SIG HT (NLOS) IN CDMA MOBILE 
Cg^UMCATION SYSTEM " — * 

Field of the Technology 

The invention generally relates to mobile communication technology, especially 
to a method that identifies whether a channel is a Line Of Sight (LOS) channel or a 
Non-Line Of Sight (NLOS) channel for a mobile station location in third generation 
cellular mobile communication system. 

Background of the Invention 

In a Code Division Multiple Access (CDMA) mobile communication system, an 
important way to raise the estimation accuracy of mobile station location is to identify 
whether the channel is a LOS channel or a NLOS channel, since the delay estimation 
error caused by a NLOS channel is the most serious one of those factors that affect the 
estimation accuracy of a mobile station location. The LOS channel can be understood 
as the transmitter is line of sight with the receiver, and the NLOS channel can be 
understood as the transmitter is non-line of sight with the receiver. 

M. P. Wylie first proposes technique that identifies a NLOS used for mobile 
station location in the article named "The Non-Line of Sight Problem in Mobile 
Location Estimation" [1]. In Sept. 9 th , 1997, he proposes a patent (with number: US 
5,974,329) to the USPTO named "Method and System for Mobile Location 
Estimation" [2], which aims at the NLOS error correction in the mobile location 
estimation and takes the NLOS identifying technique as a main element of the 
invention. In his invention the main thinking to identify a NLOS is as follow. 

(1) Make long time recorders for distance between a UE (User Equipment) and a 
base station measured in every base station; 

(2) Make smooth processing for the amount of recorded data; 

(3) Identify a NLOS with the fact that the variance measured in a NLOS (after 
smooth processing, it is caused by geomorphology feature) is much greater than the 
variance measured in a LOS (caused by system measurement error). 
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The NLOS identifying method proposed by [1] and [2] need to use the time 
correlation properties of a mobile station, i.e., the method is suitable for identifying 
LOS channel of a mobile station in a moving state, but it is unsuitable for identifying 
a LOS channel of a mobile station in a static state; furthermore, the method need to 
accumulate data for a longer period of time, which will cause a longer delay and is 
difficult to satisfy response time requirement of the FCC. 

A patent (number CN 01 105808.0) named "Method and Device for Identifying a 
NLOS channel in a CDMA cellular communication system" [3] has been applied by 
the applicant; in this patent, a LOS identifying method is proposed, which 
synthetically uses channel loss and multipath power profile to identify a LOS channel. 

Although the method in [3] has overcome the disadvantage of the method in [1] 
and [2] which only identify the LOS at the moving state of a mobile station, but in [3], 
the method has not used the power difference between paths, and, the ratio of the 
mean delay and delay spread f la in a same power delay profile to identify a LOS, 
also the selected characteristic parameters are not simple and clear enough, so the 
implementation is relatively complex. 

In the article:" Characterization of UHF multipath radio channels in factory 

buildings " [4] written by T. S. Rappaport, IEEE transaction on antenna and 

propagation, vol.37, no.8, pp.1058-1069, August, 1989, based on the indoor measured 
data, * I(T difference between LOS channel and NLOS channel at an indoor 
environment has been discussed, i.e., regularity of multipath power profile in LOS 
channel and NLOS channel at the indoor environment has been discussed; the 

conclusion is that less T I a shows that energy is concentrated on the first path and 
larger shows that energy is concentrated on the middle and trail part of the 
power delay profile. In reference [4], the discussion about r 1° is based on that the 
LOS distance is known and thereby the arriving time of LOS path can be computed; 
but in real the detected first path may not be the LOS path and also how to compute 

and use regularity of r 1° has not been discussed there. 
Summary of the Invention 
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Objective of the invention is to design a method to identify a LOS channel in a 
CDMA mobile communication system. With this method based on multipath 
searching (including de-spread, coherent accumulation and non-coherent 
accumulation), it can be identified that a channel is whether a LOS channel or a 
NLOS channel. 

One purpose of the invention is to provide a method to identify whether a 
channel is a LOS channel or NLOS channel according to the power difference 
between a direct path and a non-direct path on the same power delay profile. 

Another purpose of the invention is to provide a method to further identify 
whether a channel is a LOS channel or a NLOS channel by means of x I a difference 
between a LOS channel and a NLOS channel (wherein T is a mean delay of a 
multipath power profile, and °" is a root-mean-square (RMS) delay spread of the 
multipath power profile), based on the result of above method. 

A method for identifying whether a channel is a LOS channel or a NLOS 
channel in a mobile communication system, comprises: 

A. reading in a power delay profile; 

B. selecting a path with maximum power amplitude, i.e. a Maximum Path, from 
the power delay profile; 

C. estimating an average noise power and arrival time of a First Path and the 
Maximum Path; 

D. detecting a Local Maximum Path value within a searching window, and 
detecting whether a power ratio of the Maximum Path to the Local Maximum Path is 
greater than a threshold K\ 

E. detecting whether an arrival time difference between the First Path and the 
Maximum Path is less than a time interval T\ 

F. if the power ratio of the Maximum Path to the Local Maximum Path is greater 
than the threshold AT, and simultaneously the arrival time difference between the First 
Path and the Maximum Path is less than the time interval T, determining the channel 
being a LOS channel; otherwise determining the channel being a NLOS channel. 

After the step F, further comprises a step of detecting whether the NLOS channel 
determined by Step F is a LOS channel or a NLOS channel, by T * a difference of the 
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power delay profile between the LOS channel and the NLOS channel, which 
comprises: 

G. computing T t that is an arrival time difference between a i th detectable path 
and first detectable path; 

H. according to following formulas, computing mean delay T and root-mean- 
square delay spread a , 

» 

I. computing T ^ a that is a ratio of mean delay to root-mean-square delay spread 
of power delay profile; if is less than Delta(A ), determining the channel being a 

LOS channel, and if r I a is not less than Delta (A ), determining the channel being a 
NLOS channel; 

wherein T is a mean delay of a power delay profile and a is a root-mean-square 
time-delay spread of the power delay profile; wherein n is number of detectable paths, 
Pi is the i* path power, and i is 1 to n. 

Wherein the Delta is taken between 0.5 to 1. 

Wherein Step D further comprises, if there is no detectable path in the searching 
window of the Local Maximum Path, taking Theta (9) multiples of the average noise 
power as the Local Maximum Path power. 

Wherein the Theta (0) is taken 2. 

In Step D, the said Local Maximum Path is selected from a range within Alpha 
( a ) microsecond that delays the Maximum Path; the Alpha ( a ) width is greater 
than one chip. 
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In Step D, the threshold K in indoor environment, is set less than the threshold in 
outdoor environment, and the threshold in outdoor environment is set 10. 

Wherein Step D further comprises, dividing the threshold K into Kl and K2, 
wherein Kl > K2\ 

wherein Step F further comprises, if the power ratio of the Maximum Path to the 
Local Maximum Path is greater than AT/, determining the channel as a LOS channel; if 
the ratio is less than K2, determining that the channel is a NLOS channel; and if said 
ratio is between Kl and AT2, determining that the channel is an undetermined channel. 

The threshold Kl is taken 10 and K2 is taken 5. 

The time interval T is within three chips, and typically two chips. 

Comparing with the technical scheme in [3], the technical scheme of the 
invention uses the power difference between paths on the same power delay profile, 

and uses T I a to identify a LOS channel ( r is the mean delay and a is the RMS 
delay spread of the multipath power profile). 

A special feature of the invention is that the method for identifying a LOS 
channel can be used to a mobile station (terminal) either in a moving state or in a 
static state. Implementation of the method deploys the present signal processing 
techniques for a communication system, such as de-spread, coherent accumulation 
and non-coherent accumulation etc., so implementation is simple and compatible with 
various mobile communication systems. 

Simulation shows that the method which uses the power difference between a 
direct path and a non-direct path on same power delay profile can effectively 

differentiate a LOS channel or a NLOS channel, but also the method which uses T l a 
difference between a LOS channel and a NLOS channel and is based on the result of 
above method can effectively differentiate a LOS channel or a NLOS channel. When 
the above methods are involved in the location estimation algorithm, the mobile 
station positioning accuracy is obviously increased. 

The method for identifying a LOS channel may also be used to identify a channel 
in high speed data transmission. 
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Brief Description of the Drawings 

Fig. 1 shows diagrams of power difference between a LOS path and a NLOS 
path on a same power delay profile; wherein Fig. la is typical power delay profile of a 
LOS channel, and Fig. lb is typical power delay profile of a NLOS channel. 

Fig. 2 shows diagrams of f la difference between a LOS channel and a NLOS 
channel; wherein Fig. 2a is the curve obtained by taking theoretical arrival time of a 
LOS path as the reference point, and Fig. 2b is the curve obtained by taking arrival 
time of the measured First Path as the reference point. 

Fig. 3 shows a flowchart to determinate the LOS/NLOS channel based on power 
difference between specific paths. 

Fig. 4 shows a flowchart to identify a NLOS channel synthetically using power 
difference between paths and r 1° . 

Fig. 5 shows a system diagram of an embodiment of the invention which is used 
to identify a LOS channel. 

Embodiments of the Invention 

The invention will be described in more detail with reference to the drawings. 

Fig. 1 shows the power difference of a LOS path and a NLOS path on a same 
power delay profile. The abscissa is the relative delay with microsecond as a unit; the 
ordinate is the power with dB as a unit. In Fig. 1, the First_Path is the first path, the 
Max_Path is the maximum power amplitude path and the Second_Max_Path is the 

second power amplitude path; the Alpha ( a ) is an interval width and the 

LocaLMax_Path is the local maximum power amplitude path; the searching window 
for the Local__Max_Path is from the starting point to the ending point. 

Fig. la is a typical power delay profile of a LOS channel after having been 
processed with coherent accumulation and non-coherent accumulation; Fig. lb is a 
typical power delay profile of a NLOS channel after having been processed with 
coherent accumulation and non-coherent accumulation. 
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As shown in Fig. la, the characteristics of a LOS channel power delay profile are: 
at the beginning location or 1 to 2 chips delay of power delay profile, there is a 
strongest path (Max_Path) of which strength is Ricean-distributed; within interval 

Alpha (a ) that directly follows behind the strongest path, there are quasi-LOS paths 

caused by strong reflection of terrain, and amplitudes of them have no great difference 
with the amplitude of the strongest path but they are attenuated rapidly; after Alpha 

(a ), there are typical NLOS paths and the strength of them is Rayleigh-distributed, 

and the Local_Max_Path is searched in the searching window (from the starting point 

to the ending point). Measurement shows the Local_Max_Path amplitude (-30dB in 

Fig. la) is obviously lower than the Max_Path amplitude (OdB in Fig. la), and in 
general the Local_Max_Path amplitude is lower than the Max_Path amplitude more 
than 1 OdB (it is 30dB in Fig. la). The difference of these two amplitudes is a base to 
identify a LOS channel in this invention. 

As shown in Fig. lb, the characteristics of a NLOS channel power delay profile 
are: at the beginning location or 1 to 2 chips delay of it, there is the first path 

First_Path, then the strongest path Max_Path, and then an Alpha (a ) interval; after 

the Alpha (a ) interval they are a typical NLOS paths and the Local_Max_Path is 

searched in this searching window (from the starting point to the ending point as 
shown in Fig. lb). Measurement shows that there is no obviously difference between 
the amplitude of the Local_Max_Path and the amplitude of the Max_Path; this feature 
is a base to identify a NLOS channel in this invention. 

Fig. 2 shows the f/a difference of the power delay profile for a LOS channel 
and a NLOS channel (wherein f is a mean delay of a multipath power profile, and a 
is a root-mean-square (RMS) delay spread of the multipath power profile). In Fig. 2, 
the abscissa is a sample sequence number of the power delay spread, i.e. the power 

delay spread sequence number; and the ordinate is the ratio of T / cr 9 i.e. a mean 
delay/a RMS time-delay spread. 

Fig. 2a takes the arriving time of the LOS path as the starting point to compute 
the values of T I a curve for LOS channel and NLOS channel. It is seen that the 
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average value of T * a for the LOS channel is about 0.65 under the threshold (it is 1.0 
as shown in Fig. 2), and the average value of * 1° for the NLOS channel is about 1.5 
above the threshold; so taking r ^ (7 = 1 as a threshold, it is easier to differentiate the 
two kinds of channels. 

Fig. 2b takes the first path having been measured as the starting point to compute 

the values of r ^ <r curves for LOS channel and NLOS channel; the real-line shows 

the T I (J curve of the LOS channel, and the dot-line shows the T I G curve of the 
NLOS channel. Comparing Fig. 2a and Fig. 2b, it is shown that for the LOS channel 
the average value in Fig. 2b is decreased a little than in Fig. 2a (it is 0.65 or so), but 
for the NLOS channel the average value in Fig. 2b is obviously decreased than in Fig. 
2a (it is decreased from 1.5 to 1.0, decreasing about 0.5). 

Fig. 2b shows that, there is an ambiguous zone for T ^ a curve of the LOS 
channel and * 1° curve of the NLOS channel, but when the T lvalue is less than a 
specific value (it is 0.65 in Fig. 2b), it is determined that the channel is a LOS channel. 

This is the base of the invention to use T ^ a value for further identifying a LOS 
channel after having been identified by the amplitude difference between different 
paths. 

Fig. 3 shows a basic flowchart to identify a LOS channel or a NLOS channel by 
using power strength (amplitude) difference between a direct path and a non-direct 
path. 

In the first step, read in the power delay profile and extract its characteristic 
parameters; said power delay profile is obtained after the system has performed the 
coherent accumulation and non-coherent accumulation, and the outputs of the 
coherent accumulation and non-coherent accumulation are taken as a extracted basis 
of the characteristic parameters. 

In the second step, pick up the strongest path Max_Path; the Max_Path is the 
path with maximum power amplitude (power strength) in the whole power delay 
profile. 
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In the third step, evaluate the Average_Noise_Power, the 
Time_Of__First_Path_Arrival and the Ti me_Of_Max_Path_ Arri v al . With the present 
mature technique, the Average_Noise_Power is evaluated in two steps: first, take N 
(Af=5, for example) strongest paths away from the power delay profile, second, make 
average within the searching window for the power delay profile in which N paths 
have been taken away; wherein way to take a strongest path away is as follow: select 
the strongest path on the power delay profile, set values of samples in the range that 
covers three sample points at both side of the strongest path sample to zero, thereby 
one path power of the power delay profile having been taken away, which will be 
used for taking away the strongest path next time, is obtained; repeat N-l times of the 
above procedure to obtain the power delay profile in which the N strongest paths have 
been taken away (when the number of paths is less than N, some noise powers will be 
taken away, which does not affect application). 

In the fourth step, determine the Local_Max_Path that must be selected at the 
range that has Alpha (a ) microsecond delay from the Max_Path; if there is no 

detectable path in said range, take Theta (0) (such as 2) multiple of noise average 
power as the Local_MaxJPath value. The purpose to take said range that has Alpha 

(a ) microsecond delay from the Max_Path is to avoid the strong reflect path that 

follows the Max_Path in a LOS channel; the Alpha (a ) should be taken more than 
one chip width. 

In the fifth step, make detection according to the following two formulas: 
Max_Path/LocalJ\4ax_Path >K (1) 

Time_Of_First_Path_Amval-Time_Of_Max_Path_Arrival<T (2) 

Formula (1) detects whether a ratio of the strongest path and the local strongest 
path (Max_Path/Local_Max_Path) is greater than a threshold K. For the 
Local_MaxJPath, there are two possible situations: there is a local strongest path or 
there is not a local strongest path; in the later situation, Theta (0) multiple of the noise 
average power (0 = 2) is taken as its value. For indoor and outdoor environment, the 
threshold K is different; usually the indoor threshold is less than the outdoor threshold, 
and here K is taken 10 as the outdoor threshold. 
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Formula (2) detects whether the arriving time difference between the first path 
and the maximum path is less than a setting value T\ T is within three chips time and 
typically it is taken two chips time. 

In the sixth step, if formulas (1) and (2) are satisfied simultaneously, determine 
that it is a LOS channel; if they are not satisfied simultaneously, determine that it is a 
NLOS channel. - 

In real, the power delay profile is various; in some specific environments, the 
Max_Path/ Local_Max_Path of the LOS channel may overlap the Max_Path/ 
Local_MaxJPath of the NLOS channel. For differentiation, in this situation the 
threshold K is divided into two values Kl and K2> wherein Kl > K2, for example Kl 
is taken 10 and K2 is taken 5. 

When (Max_Path/ LocalJVIaxJPath) > Kl, it is a LOS channel; when 
(Max_Path/ Local_Max_Path) < K2 y it is a NLOS channel; and when Kl > 
(Max_Path/ Local_Max_Path) > K2, it is an undetermined channel. 

Fig.- 4 shows a flowchart for identifying a NLOS channel by using power 

difference between paths and x f a . The flowchart is divided into two parts; the first 
part uses power difference between paths to identify the LOS channels, i.e., the first 
part makes detection by means of amplitude difference between a direct path and a 
non-direct path; then the second part further identifies the NLOS channels having 
been identified by the first part to determine the LOS channel or the NLOS channel. 

The first step and second step detect whether the channel is a LOS channel or a 
NLOS channel by using amplitude difference; this is same as the procedure described 
in Fig. 3. 

During selection of the threshold AT, it is kept in mind that great missing 
detection probability is taken to ensure that less false detection probability. The 
missing detection means that parts of LOS channels are missed detecting into NLOS 
channels, and the false detection means that NLOS channels are detected into LOS 
channels falsely. Thereby the result is that there are some power delay profiles of LOS 
channel not fulfilling formula (1), and these LOS channels should be further detected 
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by the T ^ <J difference of the power delay profile between a LOS channel and a 
NLOS channel. This is what the third step to seventh step will be done in Fig. 4. 

Third step computes the delay of each path (detectable i th path) to the first 
detectable path (the First Path) with formula (3), i.e. computing the arrival time 
difference r. between the i* path and the First Path. 

r t = Time J3f Jth JPath_Arrival - Time_Of JPirst_Path_Amval ( 3 ) 

Fourth step computes the mean delay f and the delay spread a . Theoretically, 
the mean delay f and the delay spread a are defined as followings: 



r= ^U-i ! g (4) 



^JZ^lllL (5 ) 



V »^ ; (6) 
Wherein n is the number of detectable paths, r t is the delay of the i* path to 
arrival time of the LOS path (if it exists), and p t is the power of the i* path. 

In real, general speaking, it is impossible to measure a LOS path arrival time 
defined in the formulas. In order to compute * , the invention takes the real measured 
time of the First Path arrival time (Time_OfJFirst_Path_ Arrival-Time) as a starting 
point, and the obtained ^ is less than the theoretical mean delay obtained by taking a 
delay to which the direct distance between the base station and mobile station 
corresponds as the starting point, but in a LOS environment this decrement is very 
little. Nevertheless, in a NLOS environment there is a larger difference between the 7 
computed from formula (3) and the real ^ ; this can be seen by comparing Fig. 2a and 
Fig. 2b. 

Since ° is kept unchanged, in the NLOS environment the computed^/ lvalue 
is obviously less than the real value, but in the LOS environment the 
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computed T lvalue is not obviously less than the real value; the result is that the 

difference of T I a value between the LOS environment and the NLOS environment 
is decreased, and this will decrease the identification accuracy. 

Fifth step computes T I a value with formula (7) and determines what kind of 

the channel is according to an area where the computed r ^ a value located. 
Gama=r/<7 (7) 

Sixth step compares the Gama and Delta with formula (8) to identify a NLOS 
channel. 

Gama= tier < Delta ( A ) (8) 

Seventh step may take Delta (A ) = 0.6 (Delta may take a value between 0.5 to 

1), and determines a LOS channel for those that the*"/* 7 are less than 0.6 and a 
NLOS channel for those that the r I ° are not less than 0.6. 

Fig. 5 is an embodiment of the invention. In Fig. 5, the LOS identification 
module in the base station (NodeB) identifies whether a channel is a LOS channel or a 
NLOS channel for a result of de-spread, coherently and de-coherently accumulated; 
the identified result is sent up, through interface Iub, to the location estimation unit in 
the Radio Network Controller (RNC), which is called base station controller in GSM; 
the location estimation unit will use the identified result to improve the location 
estimation accuracy. 

Besides implementing the LOS identification in the NodeB, the LOS 
identification module may be put in the User Equipment and the identified result is 
sent to the NodeB and RNC through the Uu interface. 

The method for identifying a LOS channel and a NLOS channel in the invention 
may also be used to identify a channel in high speed data transmission. 
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